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We investigate how an externally imposed curvature influences lipid segregation on two-phase- 
coexistent membranes. We show that the bending-modulus contrast of the two phases and the 
curvature act together to yield a reduced effective line tension. On largely curved membranes, a 
state of multiple domains (or rafts) forms due to a mechanism analogous to that causing magnetic- 
vortex formation in type-II superconductors. We determine the criterion for such multi-domain 
state to occur; we then calculate respectively the size of the domains formed on cylindrically and 
spherically curved membranes. 
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Introduction Lipid membranes play important roles 
in maintaining cell integrity and intracellular traffick- 
ing lj. At high temperatures, the three main compo- 
nents of the membranes (saturated lipids, unsaturated 
phospholipids and cholesterols) form a homogeneous mix- 
ture; below a critical demixing temperature, the three 
components segregate into two coexistent fluid phases, 
a saturated- lipid-enriched liquid-ordered (L a ) phase 
and unsaturated-phospholipid-enriched liquid-disordered 
(Lj) phase. Small L domains formed on cell membranes 
are also referred to as rafts [2j. The L phase (or rafts) 
has a larger bending modulus than the Ld phase |5HS1- 
There have been a large body of studies ptflS] investigat- 
ing how membrane curvature influences the segregation 
of membrane lipids. All these studies are based on a phys- 
ical mechanism that the L phase because of its larger 
resistance to bending prefers flatter regions and drives 
the Ld phase to more curved regions, i.e., these studies 
have implicitly assumed what matters is the variation of 
curvature rather than the curvature itself. 

In this letter, we investigate how a uniform externally 
imposed curvature may influence lipid segregation on an 
L -Ld coexistent membrane. We show that, with the 
presence of a bending-modulus contrast, the curvature 
of the membrane induces the lipids to tilt away from 
the membrane normal and yields a reduced effective line 
tension between the lipid domains; when the curvature 
is larger than a certain critical value, lipid segregation 
leads to the formation of multiple L a /Ld domains (or 
rafts) of microscopic lengthscale rather than a complete 
separation of the two phases. We determine the cri- 
terion for the multi-domain pattern to occur; we then 
calculate respectively the size of the domains formed 
on cylindrically curved and spherically curved lipid 
membranes. We point out that although inspired by 
lipid systems the results obtained in this letter apply to 
any membranes that contain two distinct smectic phases. 

General Picture Nonzero curvatures of lipid mem- 
branes can be either induced by their intrinsic sponta- 




FIG. 1: Two types of lipid configurations of curved mem- 
branes. The dimers represent lipids, either lipids in bilayer- 
membranes or lipids in unilamellar systems. In (a), the lipids 
are perpendicular to the membrane surface; in (b), the lipids 
are parallel to each other. 

neous curvatures or imposed externally. We consider the 
latter case[17]- Examples of lipid membranes with an ex- 
ternally imposed curvature includes in vivo membranes 
attached to BAR domains of proteins [TB] and in vitro 
membranes attached to rigid substrates. Fig. [T] shows 
two extreme types of lipid configurations resulting from 
an externally imposed curvature. In Fig. [l|a), the lipids 
are perpendicular to the membrane surface and the en- 
ergy cost originates from the splay energy cost of lipid 
orientation; in Fig. [ljb), the lipids are parallel to each 
other, and the energy penalty is induced by the devia- 
tion of lipids' orientations from the membrane normals, 
or more microscopically, by the relative sliding of the 
lipids with respect to one another. 

In general, the lipid configuration of a curved mem- 
brane should be a superposition of the aforementioned 
two types. The corresponding elastic energy density thus 
includes two parts, the splay energy and the tilt energy 
penalty, 

f=\K{V-nf-\C{n-u)\ (1) 

where n and u are unit vectors representing, respectively, 
lipid orientation and membrane normal, V represents 
the two-dimensional derivative on a curved surface [19]. 
K is the splay coefficient, and C is a coefficient measur- 
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FIG. 2: Illustration of external-curvature-induced lipid- 
tilting in the phase-boundary region. The solid and hollow 
dimers represent, respectively, L and Ld lipids. The dashed 
line represents the membrane normal at the phase boundary. 



ing how strongly n is locked to u. The combination of 
K and C gives a characteristic length £ = WW/C, which 
is the penetration length of lipid tilting. For a single- 
phase membrane, given that £ is usually much smaller 
than the membrane size, most of the lipids arc perpen- 
dicular to the membrane surface [see Fig. [TJa)], and the 
bending modulus of the membrane is mostly determined 
by the splay coefficient K. A bcnding-modulus contrast 
in literature thus corresponds to a splay-coefficient con- 
trast. Note that we have ignored in Eq. a divergence 
term because including it would not change our results 
qualitatively. 

For an io-L^-coexistent membrane, the total energy 
includes the contributions from both phases. We will 
use subscripts o and d to distinguish the quantities of 
these two phases. For convenience in the expressions, we 
will also use the subscript i, with i = 1 representing the 
phase with smaller volume fraction and i = 2 represent- 
ing the one with larger volume fraction; in a case of no 
ambiguity, we will use a symbol without subscript to rep- 
resent the corresponding quantities of both phases. The 
total energy of a two-phase-coexistent membrane should 
also include a phase-boundary energy cost, which equals 
the product of the line tension t and the phase-boundary 
length. To minimize the boundary energy, the two phases 
on a flat membrane completely separate from each other, 
with each phase forming a large single domain. How- 
ever, for a curved membrane, such a single domain may 
become unstable and split into multiple domains. 

We now present the reason for the splitting of large sin- 
gle domains and determine the corresponding instability 
criterion. As illustrated in Fig. [2j when subjected to an 
externally imposed curvature, the L Q lipids tilt away from 
the membrane normals and become more parallel to each 
other so as to reduce the splay energy cost, given that K Q 
is larger than The anchoring energy, i.e., the second 
term in Eq. 0, constraints the tilting to occur within 
a narrow region of size ~ £ around the phase boundary. 
The energy change induced by such tilting can be com- 
puted. To simplify calculation, we assume that the size of 
the large single domains resulting from a complete phase 
separation is much larger than £ and that the imposed 
curvature is much smaller than l/£. Thus, in the afore- 
mentioned phase-boundary region, the splay V'-n^ can be 
approximated by 2H+(—ydji/dsi, where H is the mean 



curvature, Sj is the distance between the phase boundary 
and the lipid location in phase i (i.e., the arc length of 
the geodesic perpendicular to the phase-boundary line), 
and 7i is the angle between the membrane normal and 
the orientation of lipids of phase i. Note that the lipids 
now tilt in the plane defined by the membrane normal 
and the tangential vectors of these geodesies. The factor 
(— ) l is added into the expression of V • n, because of the 
opposite directionality of Si and s 2 . Substituting the ap- 
proximate expression of V • into Eq. and taking the 
functional derivative yields ji = 7h exp(— Si/£i), where 
7b represents the deviation angle at the phase boundary. 
The energy change resulting from the lipid tilting can 
then be obtained by minimizing over 7?, the sum of the 
energies of L and L^. The combination of this energy 
change and the line tension t produces a reduced effective 
line tension: 



t P =t 



2H 2 (K -K d ) 2 
VCjTo + VCdKd' 



(2) 



Eq. ([2]) clearly shows t e can be lowered either by increas- 
ing H or the contrast between K a and or by decreas- 
ing C and C d - 

The instability criterion of the single-domain state can 
then accordingly be obtained by simply setting t e = 0. 
We can further simplify this criterion to an approximate 
form: t/H 2 < K£. Given that K is of order 1CT 19 J 
for most lipid systems [T5] and £ can be assumed to be 
comparable to the membrane thickness and is thus of 
order lnm [SHj, the criterion is then t/H 2 < 10~ 28 J • m. 
Depending on the compositions and temperature, the line 
tension tona lipid membrane usually varies between lpN 
and O.OlpN (or even smaller) [1 113 H2] . For t ~ lpN, 
the critical curvature H c is thus of order 10 8 m _1 ; for 
t ~ O.OlpN (for example, in systems with presence of 
hybrid lipids [23]), H c is of order 10 m -1 . 

As t e becomes negative, the large single domains 
resulting from a complete phase separation become 
unstable and split into multiple small domains. In the 
following sections, we calculate the size of the lipid 
domains formed on cylindrically curved and spherically 
curved membranes, respectively. 

Lipid Segregation on a Cylindrically Curved Membrane 
For a lipid membrane subjected to a uniform cylindri- 
cal curvature, we assume there is translational symme- 
try along e 2 , the direction of the long axis of the cylin- 
der, and that the lipids align perpendicular to e z . In 
this case, a domain is a strip along e 2 , and the three- 
dimensional problem becomes a two-dimensional prob- 
lem. The splay of the lipid orientation can then be ex- 
pressed as V • rij = cos7i[l + d-fi/d8i]/R c , where R c is 
the radius of the cylinder and the #'s are the azimuthal 
angles representing the lipid location in a domain (see 
Fig. [3]). In the center of a domain, i.e., at 6i = 0, the 
lipids align along the membrane normal with 7,- = 0. 
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FIG. 3: 

surface. 



Cartoon of L\ and L2 domains on a cylindrical 



Substituting the expression of V' • rij into Eq. ([I]) and 
then differentiating the energy density with respect to 7^ 
yields 7, = Ai sinh (#ii£ c /£i)) where the j4's are the am- 
plitudes and the approximations £j -C R c and 7; < 1 
have been used. Approximately speaking, the deviation 
7i decays exponentially for domains of size much larger 
than £j and linearly for those of size smaller than £j. 
The amplitudes A\ and A2 are not independent from 
each other — they are related by the continuity condition 
7i (#1) = 72(^^2)1 where Of and #f are, respectively, the 
angular sizes of the L\ and L 2 domains (from the do- 
main centers to the domain edges). We then minimize 
the total energy (per unit length along the long axis of 
the cylinder) over the A's and obtain 



E = 



2tt 



t- 



(K, - K 2 f/Rl 



(3) 



which clearly shows a competition between the energy 
penalty of the line tension and the energy gain of the di- 
rector tilting induced by the difference between K\ and 
K 2 . The critical line tension t c , below which the com- 
pletely separated state becomes unstable, is given by the 
maximum value of the second term in the parenthesis in 
Eq. ( 3 ) , which is the same as the one we can obtain from 
Eq. 1 21 by setting t e = and H = l/(2i? c ). 

At t < t c , a stable multi-domain ground state is ex- 
pected. The domain size can be determined by numer- 
ically minimizing Eq. ([3| over 9\ (or #|) with the con- 
straint 01/02 = 0/(1 ~ y) applied, where cj> is the volume 
fraction of phase 1. The results are given in Fig.|4j which 
shows that the domain size decreases with t and quickly 
reduces to order £ D . For = 1/2 [see Fig. Qa)], because 
of the permutation symmetry of the two lipid phases, 
the curves describing how the domain size varies with 
the normalized line tension t' (= t/t c ) for the X's with a 
certain ratio v (i.e., K\ — vK 2 ) overlap with the curves 
for the K' s with a ratio of 1/v (i.e., K\ = Kz/v). For 
<f> < 1/2, this permutation symmetry is broken: the L Q 
domains formed on membranes with the L phase having 
a volume fraction 4> have a larger size than the Ld do- 
mains formed on the phase-permutated membranes, i.e., 
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FIG. 4: Relations between the Li-domain size (measured 
in unit £„, Ai = 8fR c /£ ) and the normalized line tension 
if — t/t c on a cylindrically curved membrane for (a) (j> — 0.5, 
(b) 4> = 0.4, (c) = 0.3, and (d) 4> = 0.01. The (blue) thick 
solid and dashed lines correspond to K\ = 2K2 and K\ = 
K2/2, and the (red) thin solid and dashed lines correspond 
to K\ — I.2K2 and K\ = K2/L2, respectively. We have set 
Ci = C2. 



membranes with the Ld phase having a volume fraction 
(j). Moreover, the size difference increases as <fi decreases 
[see Fig. |4ja)-(d)]. The origin of this size difference and 
its variation with (f> is given as follows. As can be seen 
from Fig. [2j with the presence of an external curvature, 
the L Q lipids provide tilt-driving forces while the Ld lipids 
resist tilt. Membranes with the Ld phase having volume 
fraction <f> possess more L Q lipids and thus tend to have 
more phase-boundary regions so that tilt can occur. A 
preference for boundary regions means a preference for 
smaller domains, and therefore leads to the aforemen- 
tioned domain-size difference and its variation with (j). 

In addition to the numerical results, we also present 
in the following analytical results for the small t limit so 
as to clearly elucidate the dependence of the domain size 
(i.e., 9{ and 0f) on the parameters. We consider two cases 
separately: (i) the volume fractions of the two phases are 
of same order, i.e., <f> ~ 1/2; and (ii) one phase has a much 
smaller volume fraction than the other, i.e., <f> <C 1/2. In 
the case of <f> ~ 1/2, when t <C t c , we have 0* <C &/-R c 
for both phases, and the deviation angles 71 and 72 both 
decay (approximately) linearly from the phase boundary 
to the domain centers. We thus approximate coth(x) in 
Eq. ([3]) as l/x + x/3. We then differentiate the energy 
given in Eq. Q with respect to Of, with the relation 
between Of and Q\ applied, and obtain 



3t xVSrl + r^/Cl - cf>) 



V R r .C> 



1 -r 



2/3 



(4) 



where C = C\ + C 2 (l — 4>)/<j> and r = K 2 /Ki. We men- 
tion the following three noteworthy points contained in 
Eq. Q: (i) a large contrast between the if's (i.e., large 
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|?' — 1|) gives a small domain; (ii) although the angular 
size 9f decreases as R c increases, the linear domain size 
(RcOf) increases with i? c ; (iii) at 4> = 1/2, the permuta- 
tion symmetry is conserved, viz, the value of 9\ does not 
change under the replacement r — > 1/r. We now turn to 
the case of <f> -C 1/2. In this case, we have 9{ <C £,i/R c 
but 6*| 3> &/Rc', 7i decays linearly and 72 decays expo- 
nentially from the phase boundary to the domain cen- 
ter. We thus replace coth.(#fii c /£i) and coth(#2-Rc/£2) 
in Eq. (|3| by €,\/{9{R c ) and coth(oo), respectively. Min- 
imizing Eq. (|3| then yields 



2t 



VW 2 



1/2 1 



1-r 



(5) 



which is inversely proportional to \r — 1| and independent 
of the radius R c (the linear domain size is thus linearly 
proportional to R c ). In addition, the angular size is now 
more sensitive to the variation of t than the angular size 
given in Eq. @ is. 

Lipid Segregation on a Spherically Curved Membrane 
We proceed to address lipid segregation on a spherically 
curved membrane. A complete phase separation now 
leads to the formation of a single (curved) disk of L\ 
lipids embedded in an L2-lipid sea. The instability crite- 
rion of the single-disk pattern an be obtained from Eq. ^ 
by setting t e = and H = l/R s - For a negative t e , a 
multi-domain pattern forms. In the case of <fi ~ 1/2, the 
shape of the formed multiple domains is complicated, and 
this complication makes an analytic calculation of the 
system energy inaccessible. We will study the (j> ~ 1/2 
case in a numerical approach in a future publication and 
consider here only the case of <f> <C 1/2. In this case, we 
have a pattern of multiple L\ disks embedded in an L 2 
sea. 

We first give the explicit form of the energy as a func- 
tion of the deviation angle 7 and the equation deter- 
mining how 7 varies from the center of a domain to the 
domain edge. We adopt a spherical-coordinate descrip- 
tion for the positions of the lipids, with the point of zero 
polar angle (i.e., 8 = 0) corresponding to the center of an 
Li domain, and we assume 7 only depends on the polar 
angle 9 but not the azimuth angle. Thus, for the total 
energy, we have 

E = 2nNjR4sm9? 



+R 



- 1 d9 sin 9 





+i? 2 



d9 sin ( 



cos7i(6>) 



^(V'-rnf+idsin 2 7l 
iX 2 (V'-n2) 2 +iC 2 sin 2 7 2 ]},(6) 
■-cot*sin 7i (6>)} > 
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where R s is the radius of the sphere, 9{ and Nj, = 
20/(1 — cos9\) are respectively the angular size and the 




FIG. 5: Relations between the Li-domain size (measured 
in unit f OJ Ai = 9{R s /^ ) and the normalized line tension 
t' = i/( c on a spherically curved membrane for K\ = 2K2 
[the (blue) thick solid line], K\ = K2/2 [(blue) thick dashed], 
K\ = l.2Ki [(red) thin solid] and Ki = K 2 /1.2 [(red) thin 
dashed] with the assumptions <f> <g 1 and Ci =Ci. 



total number of the L\ domains. The upper limit 9 U of 
the second integration equals to tt in the single-domain 
case; and its value in the multi-domain case will be ad- 
dressed later. Differentiating E with respect to 7, under 
the assumption that 7 is small, yields 



OH 

09 2 



3n „ (Ri 



+ esc 2 9) = 0. (7) 



We now proceed to determine the domain size, i.e., the 
size of the Li disks. Given the constraint <fi <C 1/2, the 
distance between the disks becomes much larger the disk 
size; and, in the L 2 sea, only a narrow annulus around 
the disks contributes to the system energy. Furthermore, 
because the domain size is small, the deviation angle is 
nonzero only at small 9's. We thus replace cot 9 and 
csc6> in Eq. ([7| by 1/9, and obtain ji oc Bli(bi9) and 
72 « BKi(6a0)> where b z = ^/(R s /&) 2 - 4 w RjZi, and 
BI ra (a:) and BK n (x) are, respectively, the modified Bcsscl 
functions of the first and second kind. To calculate the 
energy, we then follow a procedure similar to that for the 
cylindrical case and obtain 



9{ 

W(9l) := K x b x 9{ 



2{K X - K 2 f 

W(9f) 
BI (Mf) 

BiiCMf) 



K 2 b 



s BK (b 2 9l) 



(8a) 



(8b) 



Minimizing this energy yields the domain size. Fig. [5] 
shows the numerical results, which are qualitatively 
similar to the results of cylindrically curved membranes 
with small <f> [see Fig. |4]jd)]. 

Concluding Remarks We have shown that a curved 
two-phase-coexistent lipid membrane, compared with a 
fiat one, has a smaller effective line tension due to lipid 
tilting. Moreover, a large enough curvature leads to a 
negative effective line tension, and thus induces a state 
with multiple domains (or rafts) of microscopic length 
scale determined by the splay and anchoring coefficients. 
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We expect that the theoretical framework addressed in 
this paper, if applied to small intracellular vesicles, has 
some biological significance. 

We mention that the multi-domain-inducing mecha- 
nism addressed in this paper is different from that using 
elastic repulsions to stabilize multiple domains on a lipid 
vesicle [3J [23] , but analogous to the mechanism via which 
a large magnetic field induces the formation of magnetic 
vortices in type-II superconductors [35]. The tilt pene- 
tration length £ in this work is the counterpart of the 
magnetic field penetration length in superconductors. 
What corresponds to the superconducting coherence 
length is the correlation length of the density fluctuation 
of the various components of lipid membranes |22j . We 
have implicitly assumed this correlation length is smaller 
than £, viz., we have assumed a sharp boundary between 
the two coexisting phases. We will present elsewhere 
how an externally imposed curvature influences lipid 
segregation on membranes with this assumption vio- 
lated, for example, on membranes with the temperature 
close to a critical demixing point. 
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